Impact of graphene quantum capacitance on transport spectroscopy 
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We demonstrate experimentally that graphene quantum capacitance Cq can have a strong im- 
pact on transport spectroscopy through the interplay with nearby charge reservoirs. The effect is 
elucidated in a field-effect-gated epitaxial graphene device, in which interface states serve as charge 
reservoirs. The Fermi-level dependence of Cq is manifested as an unusual parabolic gate voltage 
(Vg) dependence of the carrier density, centered on the Dirac point. Consequently, in high magnetic 
fields the spectroscopy of longitudinal resistance (Rxx) vs. Vg represents the structure of the 
unequally spaced relativistic graphene Landau levels (LLs). Rxx mapping vs. Vg and B thus reveals 
the vital role of the zero-energy LL on the development of the anomalously wide z/ = 2 quantum 
Hall state. 



I. INTRODUCTION 

Graphene has been attracting great interest, much of 
which is linked to the behavior of its charge carriers that 
mimic massless Dirac fermions as a result of the relativis- 
tic band structure referred to as Dirac cones i"— Among 
various experimental techniques designed to probe Dirac 
cones, the measurement of quantum capacitance Cq^- — 
defined as Cq = e^(dnG/d£:F) with uq the charge carrier 
density and ep the Fermi level (e: elementary charge), 
is unique in that it directly probes the density of states 
(DOS) at the Fermi level. The effects of Cq usually ap- 
pear as a minor modification to the measured capaci- 
tance, expressed as Cmeas = (1/Cox + 1/Cq)~"^ (Cqx^ ca- 
pacitance of gate oxide), which thus becomes discernible 
only in the vicinity of the charge neutrality point (Dirac 
point) where 1/Cq > I/Cqx-^^ Using a very thin (10 nm) 
high dielectric gate insulator, Ponomarenko et at. were 
able to measure graphene quantum capacitance with 
high accuracy, which revealed the smearing of the Dirac 
cone by charge puddles and the broadening of relativis- 
tic graphene Landau levels (LLs) in high magnetic fields 
{B),^ These studies demonstrate that Cq indeed depends 
on ^p^'i^ unlike in conventional two-dimensional (2D) 
systems.^ However, with the standard setup for trans- 
port spectroscopy using a gate oxide a few 100 nm thick, 
the feature of graphene quantum capacitance refiecting 
the Dirac cone is not manifested, because it is overridden 
by Cox ; accordingly, the usual linear relation between uq 
and gate voltage Vg is routinely assumed as in conven- 
tional 2D systems. 

In this paper, we demonstrate that graphene quan- 
tum capacitance can have an impact on transport spec- 
troscopy over a wide carrier density range away from the 
Dirac point, through the interplay with nearby charge 
reservoirs. Here we elucidate the effect in high- mobility 
field-effect-gated epitaxial graphene, in which the role of 
charge reservoirs is played by interface states that ex- 
ist at the SiC substrate and in the gate insulator. The 
interplay between the ep-dependent graphene quantum 
capacitance and the nearly energy-independent interface 
state DOS modifies the usual relation no cx K — Ud 



(Ud: gate voltage at the Dirac point) into a nontrivial 
one £y — Vj:, [i.e., \nQ\ ex (V^ — Vq)'^] when the 

interface states dominate. In high 5, transport spec- 
troscopy vs. Vg consequently emerges as a manifestation 
of the energy spectrum of the unequally spaced relativis- 
tic graphene LLs. This additionally allows us to deduce 
the energy width of the extended state in the A/" = 
LL {N: the LL orbital index) from the width of the lon- 
gitudinal resistance {Rxx) peak. The mapping of Rxx 
vs. Vg and B provides an intuitive picture of the vital 
role of the zero-energy N = {) LL in the development of 
the anomalously wide v = 2 quantum Hall (QH) state, 
which has been commonly observed in low-density epi- 
taxial graphene. ^^^^ We explain our data coherently by 
using a device model containing interface states. Our 
results suggest that similar effects, although not as ap- 
parent, may be at work in the transport spectroscopy of 
graphene including nano- and heterostructures.^^ 



II. SAMPLE AND METHOD 



Our epitaxial graphene device was fabricated from 
graphene grown on 4H-SiC(0001) and has a top gate with 
a gate insulator made of HSQ/Si02 (120 nm/40 nm).^^i 
The device is a Hall-bar with width W = 5 /im and length 
L = 15 /im. For comparison, we also studied an exfo- 
liated graphene device of a comparable size, fabricated 
on a Si/Si02 (285 nm) substrate using the conventional 
methodi Transport measurements were performed at a 
temperature of 1.5 K using a standard lock-in technique 
with a current of less than 1 /iA. For epitaxial graphene, 
we have studied more than ten devices fabricated from 
several different SiC/graphene wafers. Here we present 
data taken from one sample for the consistency of the 
analysis, but we emphasize that similar results are ob- 
tained for all the samples studied. 
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FIG. 2. (Color online) (a) Schematic energy diagram of 
a top-gated graphene device incorporating two kinds of in- 
terface states below and above graphene, with DOS 71 and 
72. The charge density in the lower (upper) interface state 
is — eni(2). ^1(2) is the work-function difference between 
graphene and the lower (upper) interface state. The figure 
represents a case where A\ > 0, < 0, n\ < 0, 712 > 0, and 
14 < 0. (b) Equivalent circuit of (a). 



FIG. 1. (Color online) (a) Mobility vs. carrier density no 
for epitaxial and exfoliated graphene devices. (b),(c) |nG| 
vs. for (b) exfoliated graphene and (c) epitaxial graphene. 
The dotted lines indicate the relation expected from each gate 
capacitance. The dash-dotted and dashed lines indicate cal- 
culations for different interface state densities: 71 ^ 00 and 
72 = (dash-dotted), 71 = 5.0 x lO^^eV^cm"^ and 72 = 3.7 
(4.8) X 10^^ eV"^cm-2 for no > (no < 0) (dashed). 



III. EXPERIMENT AND ANALYSIS 

A. Results in zero and low magnetic fields 

Figure 1(a) compares the mobility \i [= l/nQei^a^a; x 
(L/VK)] of exfoliated and epitaxial graphene devices, 
plotted as a function of no determined from the Hall co- 
efficient in low fields (< 0.5 T). The two devices have 
comparably high mobilities, whereas a striking differ- 
ence is observed for the dependence of Incl- Exfo- 
liated graphene shows a normal Incl oc relation with 
a slope that is in accordance with the gate capacitance 
[Fig. 1(b)]. In contrast, the change in Incl with in epi- 
taxial graphene is much smaller than that expected from 
the gate capacitance, Incl = C'oxlV^ — ^ol/e [dotted line 
in Fig. 1(c)]. More intriguingly, the dependence of Incl 
on is nonlinear, being parabolic-like around the Dirac 
point. The symmetry of the feature around the Dirac 
point implies that it arises from an intrinsic property of 
graphene. 

The reduced action of on no indicates the existence 
of carrier trap states that accommodate a portion of the 
charges induced by the gate. The observed parabolic 
dependence Incl oc (Ig — Fd)^, combined with Incl oc 
e\ specific to monolayer graphene^ suggests that in the 
present case the unusual relation oc — holds, as 
opposed to the normal behavior no oc K — I^d- As we 



show below, this can be explained if the trap states have 
a large and nearly constant DOS. 

We developed a model containing interface states be- 
low and above the graphene layer, each of which has a 
constant DOS 71(2) and geometrical capacitance Cci(2) = 
^i(2)/<^i(2) with graphene, where dx{2) and £1(2) are the 
distance and effective permittivity between graphene and 
the interface states below (above) [Fig. 2(a)]. The lower 
interface states represent the dangling-bond states at the 
SiC substrate,— li^— while the upper interface states 
represent charge traps in the gate insulator ^21 By not- 
ing that the V^-induced charge redistribution between 
graphene and the interface states is influenced by both 
the electrostatic potentials associated with Cci(2) and the 
chemical potentials associated with their quantum ca- 
pacitances given by Cq = {Atlq, j AsyY^ and 



e"7i(2),^ we can describe the charge redistribution as fol- 
lows: 



^7l(2) 



, C'sl + Cs2 

no H 5 sy 



Co- 



1 + C^2/Cc2 e 



(^g-^D). (1) 



Here, Csi(2) = (1/C'ci(2) + 1/C'^i(2))~"'" is the series capac- 
itance of Cci(2) and C^i(2)- The corresponding equivalent 
circuit is shown in Fig. 2(b).^^ While Eq. (1) is a relation 
that generally holds for 2D systems, in graphene it has a 
non-trivial consequence because Ey is not proportional to 
no- That is, Cq in the equivalent circuit is a function of 
no, which leads to a nonlinear relation between no and 

We calculated no vs. by solving Eq. (1) with 
Ey = hvy/mi^ {v 1.15 X 10^ m/s^), di = d2 = 0.3 
nm, and (£1,6:2) = {eq^Seq) {eq: vacuum permittivity) 
As shown by the dashed line in Fig. 1(c), the data can 
be fitted well with physically reasonable values, e.g., 
71 = 5.0 X 10^2 and 72 = 3.7 (4.8) x 10^^ eV-^cm-^ for 
^ — > (< 0).^^ Here, we used the same 71 value as in 
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FIG. 3. (Color online) (a) Rxx mapping vs. Vg and B for epitaxial graphene. The white dashed lines, indicating the positions of 
half- filled LLs, are calculated using the interface state densities 71 = 5.0 x 10^^ eV~^cm' and 72 = 3.7 (4.8) xlO^^ eV"^ cm 
for no > (no < 0), obtained from the low- field data in Fig. 1(c). (b) Rxx mapping vs. Vg and B for exfoliated graphene. 
(c)-(e) Mapping of regions with integer (gray) and non-integer (white) LL fillings versus Vg and B, calculated using different 
interface state densities, (c) the same as in (a), (d) 71 = 5.0 X 10^2 eV-^cm-2 and 72 = 4.8 xlO^^ eV-^cm-^ and (e) 
71 = 72 = 0. The dashed lines represent half fillings of four-fold degenerate graphene LLs. (f) Vg positions of Rxx peaks in 
the data of (a) plotted vs. VNB for the = 1 and N — 2 LLs. (g) FWHM of Rxx peaks and the width F of extended states 
deduced for different / values, plotted vs. B for the N = and N = 1 LLs. / is the fraction of the extended states in the whole 
LL DOS. 



Refs.fTllandfTsl. We emphasize that the model in Ref . [Tsl, 
which assumes only interface states below graphene, can- 
not account for our results. A calculation with 72 = 
and 71 ^ 00 yields a Vg dependence much larger than 
that observed [dash-dotted line in Fig. l(c)].^^ 

Equation (1) demonstrates that the linear relation be- 
tween and Vg — Vb^ which usually does not hold for 
transport spectroscopy of graphene, does appear when 
the second term on the left-hand side dominates the first 
term. For the 71 and 72 values deduced above, we find 
this condition to be met when Ind <C 10^^ cm~^. Con- 
sequently, measurements in high B provide the spec- 
troscopy of the unique LL structure characteristic of 
Dirac electrons in graphene, as demonstrated below. 

B. Results in high magnetic fields 

Figure 3(a) shows a color-scale plot of Rxx measured as 
a function of Vg and B. Rxx peaks at \u\ > 2 separating 
adjacent QH regions appear as sets of unequally spaced 
parabolic curves centered on the Dirac point and emanat- 
ing from B = 0. [For comparison, we show the data for 
exfoHated graphene in Fig. 3(b).] The white dashed lines 



in Fig. 3(a) represent the positions at which the four- fold 
degenerate graphene LLs become half filled, i.e., v = Ai 
{i G Z)^ calculated from 71 and 72 deduced from the low- 
field data in Fig. 1(c). As the agreement between the ex- 
periment and the calculation shows, our model remains 
valid in high fields and captures the essence of the gate- 
and field- induced charge redistribution. In Fig. 3(f), the 
positions of the Rxx peaks in Vg, measured from Fd, are 
plotted for the N = 1 and N = 2 LLs with the abscissa 
scaled to VNB. The observed relation Vg — VBo: \/NB 
confirms that the mapping of Rxx peaks vs. Vg and B 
represents the structure of relativistic graphene LLs. 

As we show below, a more profound consequence of 
the relation — oc is that, when the Fermi level 
is located at one of the LLs, Vg — Vj:^ serves as a measure 
of the energy position of that LL not only at its center 
but also at its tail. By substituting in Eq. (1) with 
the energy of the graphene LL, En{B) = vV2NheB, 
where N is the orbital index of the highest occupied level, 
we are able to calculate uq and hence v = nch/eB {h: 
Planck's constant) as a function of Vg and B without 
unknown parameters. Figure 3(c) depicts the result cal- 
culated using 7i and 72 appropriate to our sample. In 
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FIG. 4. (Color online) (a) Range of fields over which the 
= 2 QH regime was observed to be extended to higher B in 



previous reports- 



■ overlaid on our Rxx data replotted vs. 



low- field carrier density Uq and B. Equation (3) calculated 
with 71 = 0, 7d, and 107d (7d = 5.0 x 10^^ eV"^cm"2) with 
72 = are shown, with the limit 71 ^ oo shown by the dashed 
line. Dash-dotted lines indicates the calculation using 71 and 
72 deduced for our sample, (b) Rxy and Rxx versus magnetic 
fields for different Vg. The dashed lines, linear extrapolations 
of Rxy from low fields, indicate the positions of = 2 expected 
from Uq. 



the diagram, regions with integer (non-integer) fillings 
are shown in gray (white) and the positions of half fill- 
ings {v = AN) by dashed lines. The Vg-B plane is filled 
mostly with integer- regions, and non-integer u appears 
only in narrow regions around half fillings. This happens 
because the Fermi level can lie in the gap between LLs by 
virtue of the interface states, whose DOS exists in the LL 
gap. The impact of the interface states is made clearer 
by comparing Fig. 3(c) with the calculations performed 
with smaller 72 [Figs. 3(d) and 3(e)]. As 72 increases, 
the LL fan diagram becomes more nonlinear and, at the 
same time, integer regions become increasingly dominant 
over non-integer regions. 

Our model accounts for the main feature of the exper- 
imental results sufficiently well. Yet, it is noticeable that 
at low B the measured Rxx peaks are broader than the 
calculated widths of the non-integer regions. This is be- 
cause the model does not take into account the LL broad- 
ening caused by disorder. This, in turn, allows us to ex- 
tract the energy width F of the extended states from the 
widths of the Rxx peaks. Note that F is otherwise not ac- 
cessible through other conventional techniques including 
capacitance,^ infrared, and scanning tunneling^^ spec- 
troscopy, which do not distinguish extended and localized 
states. When the overlap with adjacent LLs is negligi- 



ble, the full width at half maximum (FWHM) AVg of an 
Rxx peak is shown to be related to F as 



F = 



Co: 



Csl + Cs2 1 + Cj2/ Cc2 



I^AV,-fa-^y (2) 



Here, / is the fraction of the extended state with respect 
to the whole LL DOS and a is the fraction of the LL DOS 
within its FWHM {a = 0.761 for Gaussian and a = 0.5 
for Lorentzian).^^ In Fig. 3(g), we plot the measured AV^ 
and F deduced for the N = and N = 1 LLs assuming 
different / values as a function of 5.^^ Our analysis yields 
F ^ 16 ± 2 meV for the N = LL with reasonable 
accuracy at low B, even without any knowledge of the 
exact f or a values. Although F for the N = 1 LL 
is contingent on the stronger overlap with the N = 2 
LL, we can say that F's for the = LL and A" = 1 
LL at low B are of similar magnitude. This allows for 
an interesting contrast to be made with our results and 
the activation gap measurements reported for exfoliated 
graphene, the latter suggesting that the A' = LL is 
much sharper than higher LLs.^^ We may also make a 
nice comparison with theory predicting that the A" = 
and other LLs exhibit different behavior for certain types 
of disorder. ^^^^ 

As shown in Fig. 4(b), the u = 2 QH effect persists up 
to unexpectedly high fields, a feature commonly observed 
in epitaxial graphene with low carrier densities ^i^Tii^ The 
effect has recently been discussed by Janssen et al. in 
terms of 5-induced charge transfer from the dangling- 
bond states to graphene. As we have demonstrated in 
Figs. 3(a) and 3(c), when the interface state DOS is high, 
Rxx mapping vs. and B mimics the graphene LL struc- 
ture. This intuitively shows that the anomalously wide 
v = 2 QH region is a consequence of the special property 
of the A^ = LL, namely, that it has zero energy inde- 
pendent of B. With our model, the critical field at which 
the Fermi level enters the A" = LL is given by 



Be 



h 
2e 



(Csl + Cs2 



(3) 



as a function of the zero-field carrier density Uq?^ The 
second term represents the effect of 5-induced charge 
transfer. Equation (3) reveals that interface states both 
above and below the graphene layer contribute to the B- 
induced charge transfer, which increases Be irrespective 
of whether they are donors or acceptors. Equation (3) 
also allows us to analyze all reported data on the anoma- 
lously wide u = 2 QH state, including ours, in a single 
framework by plotting them vs. Uq. 

Figure 4(a) summarizes the range over which the 
ly = 2 plateau was extended toward high B in previ- 
ous reports, ^^^^ overlaid on our data replotted vs. tiq. 
The figure also shows Eq. (3) calculated with 72 = and 
different 71 values, with the Hmit 71 ^ 00 shown by the 
dashed line. The calculation with 7's for our sample is 
shown by the dash-dotted lines. It must be noted that 
the onset of transport in the A^ = LL was not observed 
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in any of the previous reports; the high- 5 end shown in 
Fig. 4(a) only indicates the maximum field used in the 
measurements. These fields, which fall within the range 
from 7i 7d to 107d (7d = 5 x 10^^ eV~^cm~^) if 72 = 
is assumed, thus gives only a lower bound for 71. More 
importantly, as Eq. (3) shows, any type of doping that 
modifies the carrier density in graphene at 5 = can 
result in 5-induced charge transfer in high B and thus 
increase Be. It is therefore possible that molecules or a 
gate insulator used to reduce the B = density in these 
studies^^^^'^^^^ is the main cause of the anomalously 
wide v = 2 region. Since the contributions of 71 and 72 
are not distinguishable from 5-sweep measurements as 
evident from Eq. (3), it is essential to examine Vg de- 
pendence, as we have demonstrated here, to identify the 
exact origin of the anomalously wide u = 2 QH region 
and to predict the onset of transport in the = LL. 



IV. CONCLUSION 

In summary, we have demonstrated that the graphene 
quantum capacitance can have an impact on transport 
spectroscopy through the interplay with the interface 
state DOS. With high interface state densities, gate- 
sweep transport measurements serve as energy spec- 
troscopy of graphene. Our model and analysis will be 
useful for understanding various transport spectroscopy 
results for gated graphene devices, not limited to epitax- 
ial graphene but also including exfoliated graphene and 
their nano- and heterostructures^^ 

We thank M. Ueki for help with device fabrication. 
This work received support from Grants-in-Aid for Sci- 
entific Research (Nos. 21246006) from the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan. 



Appendix: Derivation of the model 

1. Charge equilibrium in a gated graphene device 
containing interface states 

We first derive an equation describing the charge equi- 
librium in a gated graphene device in the presence of in- 
terface states both above and below the graphene that act 
as charge reservoirs. For epitaxial graphene, the interface 
states above and below the graphene correspond to the 
interface states in the gate insulator and the dangling- 
bond states at the SiC, respectively. We focus on a case 
where the charge carriers are electrons. The model can 
be easily generalized to a case where the charge carriers 
are holes. 

Our model assumes two kinds of interface states with 
a constant density of states 71 and 72 located below and 
above graphene at distances of di and ^2, respectively. 
[Hereafter we use subscript 1 (2) to specify quantities as- 
sociated with the lower (upper) interface states.] These 



interface states are capacitively coupled with graphene 
via the geometrical capacitance Cci(2) = ^i(2)/<^i(2)7 
where e^2) is the effective permittivity between the lower 
(upper) interface and the graphene, and thus act as 
charge reservoirs. Due to the work- function differences 
Ai and A2 between graphene and the interface states, 
charges are redistributed among them until their electro- 
chemical potentials equilibrate, where the interface states 
act as donors or acceptors depending on the signs of Ai 
and 742. Gating, which can be either via the field ef- 
fect or molecular /polymer doping, induces the charges 
to redistribute while maintaining the equilibrium among 
the electrochemical potentials of the interface states and 
graphene. 

We express the equilibrium charge densities in the 
lower interface states, graphene, upper interface states, 
and gate as — eni, — enc, — en2, and — engate, respec- 
tively, where e (> 0) is the elementary charge. Charge 
conservation requires that 

^1 + + n2 + ngate = (Al) 

The equilibrium between the electrochemical potentials 
of the lower interface states and graphene is expressed as 

Ui + ^ +^ni=£F, A2 
V 71 / ^1 

where Sy is the Fermi level of graphene. The second 
term on the left-hand side represents the change in the 
chemical potential of the interface states caused by the 
addition or removal of electrons. (That is, the sum of 
the two terms in the parenthesis represents the chemical 
potential of the lower interface states.) The third term 
represents the electrostatic potential difference between 
the lower interface states and graphene resulting from the 
transfer of charge —eni through the plane between them 
(i.e.. Gauss's law.) A similar relation holds for the upper 
interface states and graphene 

M2 + — + (n2 + ngate) = s:f, (A3) 

V 72 / ^2 

where the terms describing the electrostatic potential dif- 
ference between the upper interface states and graphene 
result from the transfer of the charge — e(n2 + ^gate) 
through the plane between them. 

By introducing the quantum capacitances ^^1(2) = 
e^7i(2) of the interface states along with the geometrical 
capacitances Cci(2) = ^i(2)/<^i(2)7 we are able to deal with 
the effects of the chemical potentials and the electrostatic 
potentials on an equal footing and rewrite Eqs. (A2) and 
(A3) in the following more intuitive forms 

M + (7^ + 7^ ) e'ni = £F, (A4) 

\Ocl ^71/ 

. M M 2 

A2 + 7; V e 712 + T^^gate = ^F- (A5) 
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By eliminating rii and 722 from Eqs. (Al), (A4), and 
(A5), we obtain the relation between ngate and uq 



HQ 



^gate 



1 + C^2/Cc2 e2 

(A6) 

where Csi(2) = (1/^^01(2) + 1/C'^i(2)) is the series ca- 
pacitance of Cci(2) and C^i(2). By substituting = 
hv^TTUG {h: Planck's constant divided by 27r, v: the 
Fermi velocity) into Eq. (A6), we obtain uq as a func- 
tion of ngate- 



2. Gate-voltage dependence of carrier density 



With field-effect gating, ngate is related to the voltage 
applied to the gate. As explained later, the approxi- 
mate relation ngate ^ —Co^Vg/e holds, where Cox is the 
capacitance of the gate insulator. Thus, by substitut- 
ing ngate ^ —Co^Vg/e together with = hv^/nriG into 
Eq. (A6), we obtain uq as a function of Vg. 

By setting uq = and = in Eq. (A6), we find 
that the gate voltage Vd at the Dirac point can be given 
in the form 



■ CS2A2 



eCo: 



(A7) 



Hence, the relation between Vg and uq is expressed as 
follows, 



Co: 



(^g - ^d) 



Cj2 
Cs2 



nG + (Csl+Cs2) 



hv^TTUQ 



(A8) 

It should be noted that the parabolic dependence Incl cx 
(l^ — Fd)^ occurs when the second term on the right-hand 
side dominates the first term. Note also that the effects of 
the upper and lower interface states are not symmetric 
in the equation. As explained below, this refiects the 
location of the interface states with respect to graphene 
and the gate. 



3. Equivalent-circuit model 

The interplay between the interface states and 
graphene can be represented in a more transparent and 
intuitive manner with the help of an equivalent-circuit 
model. The equivalent-circuit model can be constructed 
by considering the differential relations among the vari- 
ations of the charge densities in the interface states and 
graphene, — edni, — edn2, and — ednc, induced by a 
small change in the gate charge — edngate- By taking the 
derivatives of Eqs. (A4) and (A5) and using the relation 



1/Cq = (l/e^)d£:F/d^G, we have 



1 




1 




edni = 










1 




1 




edn2 = 




C^2 







edriG ( = Vi), 



(A9) 



ednc 



1 



c2 



■e (dni + dnc) ( = V2). 



(AlO) 



When deriving Eq. (AlO), we used the charge conserva- 
tion [Eq. (Al)]. Eq. (A9) implies that capacitors with ca- 
pacitances Csi and Cq connected in parallel hold charges 
edni and ednc, respectively, at a common bias voltage 
Vi. Eq. (AlO) indicates that the capacitor with C^2 
holding charge edn2 is connected in parallel to the se- 
ries capacitors with Cq and 0^2 holding charges ednc 
and e (dni + dnc), respectively, while sharing a common 
bias voltage V2. With these considerations, we are able 
to deduce the equivalent circuit as shown in Fig. 5. 



V2- 



V,- 

Cs1 



Cc2 

HI- 



-e{dn^+dnQ) +e{dn^+dnQ) 



-ednQ +e6nQ 



Cy2 



+edngate -ed/igate 



-edn^ +edn9 



FIG. 5. (Color online) Equivalent circuit for the differential 
relations in Eqs. (A9) and (AlO). The model describes the 
changes in the amount of charge stored in each capacitor in- 
duced by a small change in the gate voltage. 

The equivalent-circuit model also tells us that the to- 
tal gate capacitance Ctot, which relates the applied gate 
voltage and the induced gate charge via dngate/dlg = 
-C'tot/e, is given by 



1 



Cu 



1 

Coy 



1 



(All) 



In the absence of interface states, the above equation 
reduces to Ctot~^ = which is the relation that 

applies to usual cases. On the other hand, if the interface 
state densities are very high such that Cox <C C^i, C^2 
as demonstrated in our experiments, Eq. (All) reduces 

to Ctot ^ C'ox- 
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4. Landau-level filling in high magnetic fields 

In this section, we derive equations that represent the 
way in which the Landau-level (LL) filling in graphene 
varies as a function of Vg and magnetic field B in the 
presence of interface states. As shown in the previous 
section, in zero field, uq and Vg follow the relation, 



Co: 



(^g - ^d) 



C^2 
Cs2 



In high fields, the energy spectrum of graphene is quan- 
tized into di screte L Ls whose energies are given by 
En{B) = vV2NheB, where (= 0, 1, 2. . . ) is the 
LL orbital index. When the Fermi level is located at one 
of the LLs, the Fermi energy is just the energy of the 
highest occupied LL, and thus in Eq. (A12) is given 
as a function of B SiS = En{B) = vVWheB, with 
N being the orbital index of the highest occupied LL. 
Noting the relation uq = veB/h {h: Planck's constant) 
and inserting and uq into Eq. (A12), we obtain the 
relation for K, i/, and N as follows. 
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vVWheB 



(A13) 

When the Fermi level resides in the A'th LL, the u 
value spans a range from 4(A^ - 1/2) to 4(A^ + 1/2) as B 
or Vg is varied, where the factor 4 arises from the spin and 
valley degeneracy of graphene LLs and the term 1/2 from 
Berry's phase. Thus, by substituting z/ = 4(A/' — 1/2) 
or 4:{N -h 1/2) into Eq. (A13) and solving it to obtain 
B for each value of V^, we are able to map the region 
where the A'th LL is partially filled. A result of such 
calculations is shown in Fig. 6. The white regions are 
those where the highest occupied LL is partially filled, 
i.e., where u is non-integer. As the figure demonstrates, 
these regions with non-integer u are separated by broad 
regions corresponding to integer z/ (shown in gray). As 
discussed in the main text, these regions emerge as a 
consequence of the large interface state densities. 

The dashed line along the center of each white region 
represents the point where the Nth. LL becomes exactly 
half filled, and this is obtained by solving Eq. (A13) with 
1/ = 47V. Note that these dashed lines, which would cor- 
respond to peaks in Rxx, are unequally spaced in V^, 
with the distance increasing with decreasing lA"]. Con- 
sequently, at each B, the u = 2 regions are much wider 
in Vg than any other integer region. This results from 
the fact that graphene LLs are unequally spaced. When 
the interface state densities are very high, the second 
term in Eq. (A13) dominates the first term, leading to 
|Vg -Vb\ oc VNB for |A^| > 1. As a result, the LL fan 
diagram observed in the mapping of Rxx on the Vg-B 
plane mimics the energy diagram of graphene LLs plot- 
ted versus B. 

A noteworthy exception to the above argument occurs 
for the A" = LL. Because of its zero energy, the sec- 
ond term in Eq. (A13) vanishes, leaving a linear relation 




FIG. 6. Mapping of regions with integer (gray) and non- 
integer (white) Landau level fillings as a function of Vg and 
calculated with 71 = 5.0 x 10^^ eV"^cm"^ and 72 = 3.7 x 10^^ 
eV~^cm~^. The dashed lines represent half fillings of four-fold 
degenerate graphene Landau levels. 



between Vg and 5. As a result, the line defining the high- 
field end of the z/ = 2 region appears as a straight line on 
the Vg-B plane, in marked contrast to other lines separat- 
ing integer and non-integer regions that are all curved. 
As detailed below, this unique property of the A" = 
LL is responsible for the anomalously wide u = 2 region 
often seen in the magnetic-field-sweep data of epitaxial 
graphene. 



5. Comparison with experiments with fixed carrier 
densities and onset of transport in the N = LL 

The above discussion has focused on how the LL fill- 
ing in graphene varies as a function of and B in the 
presence of interface states. However, in many cases, ex- 
periments on epitaxial graphene in high magnetic fields 
employ samples with static gating, such as polymer or 
molecular gating. The above formulation derived as a 
function of Vg is not directly applicable to these ex- 
periments. To allow experiments with different gating 
schemes to be compared in the same framework, in the 
following, we derive a formula that does not contain Vg 
and, instead, contains the carrier density in zero field 
{uq) as the input parameter. 

As mentioned above, Ctot can be regarded as a con- 
stant that does not depend on B. Accordingly, for a 
given value of V^, the left-hand side of Eqs. (A12) and 
(A13) can be regarded as a constant independent of B. 
This allows us to eliminate Vg by subtracting Eq. (A12) 
for 5 = from Eq. (A13) for a finite B. We obtain the 
following equation. 



Csl + Cs2 



ueB 



ANeB 



(A14) 

where we used the relation = hv^/wn^ that holds at 
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B = 0. When the interface state density is very high, 
neglecting the first two terms in Eq. (A14) yields B = 
hn^/ANe. This equation indicates that, even when the 
interface state density is very high, the magnetic field 
positions of Rxx peaks v = AN {N > 1) follow exactly 
what is expected from the zero-field carrier density. 

The above argument is not applicable to the N = LL. 
The center of the Rxx peak associated with the N = 
LL is fixed at the charge neutrality point and therefore is 
not a function of B. We thus take the critical field B^ at 
which the Fermi level enters the N = LL as a measure 
characterizing the width of the u = 2 region along the B 
axis. By setting u = 2 and N = in Eq. (A14), Be is 
obtained as a function of \ 



Be 



h 
2e 

h 
2e 



(Csl + Cs2, 



(A15) 



n^ + (Csi+Cs2)' 



-0 1 



where = hvyirn^ is the Fermi energy of graphene at 
5 = 0. Note that the first term of Eq. (A15), hn^^/2e 
(=5^), is simply the magnetic field at which v = 2\s ex- 
pected from the zero-field density. The second term rep- 
resents the effects of interface states that transfer charged 
carriers to graphene as B is increased, thereby pinning 
the filling factor at = 2 over an unexpectedly broad 
range of B. In Eq. (A15), the second term can be viewed 
as the density of carriers stored in reservoirs with the 
capacitances C^i and Cs2, which are both biased at a 
voltage Sp/e at B = 0. The reservoirs here are the in- 
terface states, whose capacitances are determined as se- 
ries capacitances of their quantum capacitances and the 
geometrical capacitances with respect to graphene. In 
high magnetic fields, all of these charges are transferred 
to graphene when the Fermi level enters the zero-energy 
N = {) LL, where the "bias voltage" vanishes. 

It is also worth noting that, when expressed as a func- 
tion of tiq as in Eq. (A15), the effects of the upper and 
lower interface states on B^ appear to be equivalent, and 
are only weighted by C^i and 0^2 • This clearly shows that 
the individual contributions of the upper and lower inter- 
face states cannot be distinguished from a 5-dependent 
study for a fixed value of tlq. This contrasts with the 
case when the equations are expressed as a function of V^, 
where the effects of the upper and lower interface states 



appear asymmetric, reflecting their asymmetric locations 
with respect to the gate and graphene. This allows us for 
the first time to estimate the individual contributions of 
the upper and lower interface states, which highlights the 
significance of our 1^-dependent study. 



6. Energy width of extended state 

Our model has so far not taken into account the broad- 
ening of LLs due to disorder. When the overlap with 
adjacent LLs is negligible, effects of finite LL widths can 
be incorporated in the model as follows. Suppose that 
the gate voltage is swept from VgioV^^ 5Vg so that the 
Fermi level moves from to £y + within a disorder- 
broadened LL across its center. The resultant change in 
the carrier density from uq to uq -|- Suq is related to 5Vg 
and 3E through Eq. (A12) as follows 



H ^ 



Co 



1 + C^2/Cc2 e 



■8V^. (A16) 



If 5E corresponds to the full width at half maximum 
(FWHM) of the disorder-broadened LL, Stiq can be ex- 
pressed as 5nQ = a-AeB/h^ where a is the fraction of the 
LL density of states within its FWHM [a = erf (In 2) ^ 
0.761 for Gaussian and a = 1/2 for Lorentzian] and the 
factor 4 originates from the spin and valley degeneracy. 
We thus obtain 



4eB Oi+C^ 

a ■ — 1 ^ oE 

h 



Co 



1 + C^2/Cc2 e 



■W^. (A17) 



Since transport measurements probe only extended 
states, we introduce a parameter / that denotes the frac- 
tion of the extended states in the total LL density of 
states. The width Al^of the longitudinal resistance peak 
can thus be related to the energy width F of the extended 
states as 



Co: 



C'sl + Cs2 1 + Ce2 



AV,-fa-^) (A18) 



by replacing a, 5E^ and SVg in Eq. (A17) with /a, F, 
and AV^, respectively. 
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Within our model, Vb contains the information on the 
work-function difference ^1(2) between graphene and the 
lower (upper) interface states (see Section 2 of Appendix.). 
We can thus deduce Ai = - 1.0 - 1.3 eV, A2 = -0.08 - 
—0.13 eV for Vb = —33 V, by applying our model to the 
situation before and after the gate insulator is deposited. 
Note that the electron density in our epitaxial graphene is 
reduced significantly from ^ (3 — 4) x lO^^cm"^ to several 
lO^^cm"^ simply by the deposition of an HSQ/Si02 insula- 
tor. This is similar to the methods used to reduce electron 
densities^*^"^^ or even to improve sample quality.^^ 



